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We consider prescriptions that are free from the direct charge-screening effects by quark loops 
and enable us to clarify the confining nature of a vacuum. We test two candidates for an order 
parameter, a Polyakov loop (P) evaluated in zero-triality backgrounds and fermionic determinants 
(©1,2) with non-zero triality. Especially, 23i,2 has very small fluctuations in comparison with a 
Polyakov loop in zero-triality sector, and seems to well reflect the characteristic of a vacuum. Such 
prescriptions could be still usable for the clarification of the confinement property of a vacuum. 
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I. INTRODUCTION 

Clarification of the chiral phase transition and the con- 
finement/deconfinement transition are the longstanding 
central issues in QCD. There have been numerous efforts 
to understand the phase structure of QCD at zero and 
finite temperature or density regions []|. Chiral phase 
transition is detected by monitoring the order param- 
eter, chiral condensate {4'fp), whereas a Polyakov loop 
is widely used in order to know whether a system is in 
the confinement phase or not. The expectation value of a 
Polyakov loop is zero when a system is in the confinement 
phase, and it can be finite in the deconfinement phase. A 
Polyakov loop may be considered as an order parameter 
for the color confinement. These two "order parame- 
ters" unfortunately do not go together. (Interestingly, a 
Polyakov loop and Dirac eigenvalues, which also refiect 
the chiral phase transition, can be related [3, Q-) While 
chiral condensate ('0V'} can be a good order parameter 
if we neglect the small current quark masses, a Polyakov 
loop can serve as an order parameter only in pure gauge 
systems without quarks (or with infinitely heavy quarks), 
where deconfinement transition is expressed by the spon- 
taneous Za-symmetry breaking. The presence of dy- 
namical quarks explicitly breaks the Z3 symmetry and 
a Polyakov loop results in an approximate order parame- 
ter. The color fields induced by a Polyakov loop is readily 
screened by the dynamical quark loops that twist along 
the Euclidean time direction. Nevertheless, a Polyakov 
loop surely reflects the confinement nature and has been 
used in many studies. 

It may be however desired to clarify the confinement 
nature in a clearer manner. For example, the chiral sym- 
metry breaking surely has an impact on the properties of 
quarks. Current quarks would change their natures due 
to the spontaneous chiral symmetry breaking acquiring 
the large effective masses [J, Q , and therefore a Polyakov 
loop screened by "constituent" quarks could be affected 
by the chiral phase transition. The apparent coincidence 
of the critical temperatures for the chiral and confine- 
ment transitions is still under debate. Even if we eval- 
uate a Polyakov loop correlator (P(O)P(x)), it is not so 
straightforward to clarify the color-confinement nature 
since we always have string-breaking effects. In high- 



density systems, the concept of color confinement could 
be obscure. Quarks may freely move from one hadron 
to another in sufficiently dense hadronic systems, even if 
the vacuum is still in "color confinement" phase. In fact, 
"quarkyonic" phase has been proposed recently p, 0], 
where quarks are confined but their degrees of freedom 
can dominate the system. In any case, as long as we insist 
on a Polyakov loop, it is needed to single out the vacuum 
property independently of the direct charge-screening ef- 
fects. Search for order parameters have a long-standing 
history, and several studies have been performed so far 
[E B [l3, \1M ■ We consider prescriptions that remove the 
dynamical quark loops which directly screen the color 
fields from a Polyakov loop. 



II. FERMIONIC DETERMINANT 

We introduce the basic ideas [3, M, [ifl \lM i'^ this 
section to make this paper self-contained. We assume 
SU(3) lattice gauge theory with single quark-flavor in 
what follows. (Extensions to other cases are simple.) The 
QCD partition function is expressed as 



Z= I dUdetD e^^sl^l, 



(1) 



with D a (lattice) Dirac operator, in which non-zero 
chemical potential can be introduced, and Sg[U] being 
the gauge action. In advance of the gauge fields' integra- 
tion, quark fields are integrated out, and all the quark 
dynamics is encoded in det D. The quark loops that can 
directly screen the color fields from a Polyakov loop orig- 
inate from such a fermionic determinant. 

This fermionic determinant can be expanded in terms 
of two types of quark loops (See Fig. [1]); ordinary quark 
loops and wraparound quark loops. Ordinary loops ex- 
ist entirely in a system and can be smoothly shrunk. 
Wraparound loops twist along the imaginary-time direc- 
tion, and they are the very loops that can directly screen 
a Polyakov loop. Such wraparound loops are also respon- 
sible for the finiteness of quark density at finite chemical 
potential, and are important ingredients. 

A fermionic determinant can be generally divided into 



detD = 




wraparound loop 



FIG. 1: Schematic figure for the expansion of a fermionic 
determinant in terms of quark loops. 



three terms; 



V = dctD= Y^ 



V, 



(2) 



with Vi {i — 0,1,2) the terms which contain (3fc + i) 
wraparound quark loops, respectively. Here k is an 
integer and can be negative, and a wraparound anti- 
quark loop is counted as —1. One connected quark loop 
which twists m-times along the Euclidean time direc- 
tion is counted as m quark loops. The total number 
of quark loops is defined as a net number: In case a con- 
tribution contains as many wraparound quark loops as 
wraparound anti-quark loops, the number is defined as 
zero. We concentrate only on wraparound loops, since 
ordinary quark loops are irrelevant in the following ar- 
gument. Next, we consider the uniform Z^ transforma- 
tion; Ut{t = 0) -^ zUt{t = 0) with z = exp(2/37ri). 
Here Ut{t = 0) are the link variables on a lattice. By 
such a transformation, a Polyakov loop P is rotated as 
P -^ zP. This Z3 transformation also affects dynami- 
cal (wraparound) quark loops, and "Di are transformed 
as Vi -^ z^T>i. Using this property, we can single out 
each Vi by means of such Z3 link-variable transforma- 
tions. Denoting the fermionic determinant after n-times 
such Z3 transformations as P(7i), Vi can be obtained as 



V, 



n=0,l,2 



V{n) X z 



(3) 



(See earlier works [^, M, U^ UM ) With vanishing quark 
chemical potential, Vq = Vq and V^ = V2 hold due to 
the charge conjugation symmetry. As we have mentioned 
above, Pj contains (3fc-f i) wraparound quark loops, and 
hence what screens the Polyakov loop's color fields are 
the quarks encoded in I?2- In fact, the complex phase 
associated with V2 rotates a Polyakov loop into the real 
sector. (2?o f^nd I?i do not.) The important note is that 
quarks in Vq and Vi cannot screen a fundamental charge 
completely. Then, if we evaluate only Vq, the response 
of a Polyakov loop would be the same as that obtained 
in a quenched system. Vq is in fact Za-transformation 
invariant and does not cause the explicit Za-symmetry 
breaking. We note here that this removal of the ex- 
plicit Za-symmetry breaking is essentially different from 
the quark-loop quenching. Evaluating Vi corresponds to 
computing the partition function Zi in the i-triality sec- 
tor (T = i), 



Z^^ 



E 

|3fc+j> 



{■ik + i\e-'^" \'ik + i). 



(4) 



Here, |3fc -I- i) denote all the possible states that con- 
tain (3fc -I- i) net quarks {k is an arbitrary integer). The 
vacuum still contains dynamical quark loops, and if the 
quark chemical potential is finite, there will exist ik net 
quarks (fc net "baryons") in the system. String-breaking 
phenomena caused by dynamical quarks will be also re- 
produced. It was also shown in Ref. [il| that such a 
projection does not change thermal properties. We ex- 
pect that it would be possible to unambiguously clarify 
the confining property of a vacuum independently of the 
direct charge-screening effect by dynamical quarks. (If 
we employ an anti Polyakov loop, the roles of Vi and V2 
switch positions with each other.) 



III. NUMERICAL TESTS 
A. Polyakov loop 

In this subsection, we individually and explicitly in- 
vestigate the distributions of a Polyakov loop evaluated 
with Vi{i = 0,1,2). We generate quenched gauge con- 
figurations with the plaquette action at /? = 5.7 on 4** 
lattice, and perform reweighting using Vi as well as the 
full determinant V generated with the Wilson fermion at 
K = 0.1600. 

Fig. [2IJ upper) and Fig. [2lJlower) are the Polyakov-loop 
distributions in the complex plain evaluated with V and 
Vq, respectively, B.t (3 — 5.7, where the broken Z3 sym- 
metry is observed. As expected, the strength in the real 
sector is much enhanced by the effect of V, which can be 
seen in Fig. [21 upper). On the other hand. Fig. [^Jlower) 
shows the similar strength to the quenched case, which 
was expected in the previous section. Especially one can 
find a three-peak structure in the complex plain. 

We show the same plots obtained at /3 = 5.0 in Fig. [31 
where no Z3 breaking can be found. The strength in the 
real sector is again (but weakly) enhanced by the effect of 
V as can be seen in Fig. [3]^ upper), whereas Fig. [3](lower) 
shows the similar distributions to the quenched system. 
This prescription, where we measure a Polyakov loop in 
zero-triality backgrounds (evaluating Vq), is free from 
the direct screening effects or the string-breaking effects 
by quark loops, and seems to enable us to clarify the 
confining nature of a vacuum. 

For further clarification, we present the distributions 
reweighted with Pi_2 in Fig-Hl One can observe that only 
V2 can rotate a Polyakov loop into the real sector, which 
implies that directly-screening quarks are coming from 
I?2- The Polyakov loops evaluated with Vi are rather 
uniformly distributed around the origin. Then, Vi or 
Vq 4- Vi may be also employed for this prescription. The 
nonvanishing Polyakov loops are sometimes related to the 
explicit Za-symmetry breaking caused by fermions. In- 
cluding Vi actually breaks the Z'3 symmetry explicitly, 
but it gives a similar Polyakov loop distribution to the 
Z3-symmetric quenched case, which is natural if we take 
into account that quarks in Vi cannot screen a Polyakov 







FIG. 2: The distributions of the Polyakov loops evaluated 
with ©(upper) and X'o(lower) at /3 = 5.7, where the broken 
Zz symmetry is observed. 



FIG. 3: The distributions of the Polyakov loops evaluated 
with ©(upper) and ©o (lower) at /3 = 5.0, where no Z-j, sym- 
metry breaking is observed. 



loop completely and hence the situation is physically sim- 
ilar to the quenched system. 

So far, the response of the Polyakov loop evaluated 
with Po has been simply similar to the quenched case, 
and the difference is unclear. In order to show the differ- 
ence, we perform a test with non-vanishing quark chem- 
ical potentials. We hereafter focus on the phase angles 
of Do and T) at finite chemical potential. (Averaging the 
fermionic determinants at finite chemical potential was 
suggested and performed some years ago in the context 
of the reduction of phase fluctuations |12. flSl. \\M^ . ) 

Fig. [5] shows the phase angles of Vq and T> obtained at 
/3 = 5.7, K = 0.1640, and < /i < 0.5, with one gauge 
configuration. The phase angle of the full determinant 
T) grows as we increase the chemical potential /i, which 
is caused by the asymmetry between quarks and anti- 
quarks. However, as can be found in Fig.[5l the phase an- 
gle of I?o remains zero in the small-/i region, and rapidly 
grows above some value of \x. The reason for this behav- 
ior would be that the number of wraparound quarks in 
■Do is 3fc. When /i is small enough, it is not likely for the 
system to accommodate three net quarks, and Do is still 
symmetric between quarks and anti-quarks, which im- 



plies fc ^ 0. The phase angle at small- /j, region is brought 
about by Di and D2. At large /i, the charge conjugation 
symmetry in Do is largely broken, and the phase angle 
can be finite, which implies /c > 0. 



B. Distributions of Fermionic determinants 

We next take a look of the distributions of fermionic 
determinants themselves in quenched QCD in this sub- 
section. Fig. [6] shows the distributions of T)\ in the 
complex plain computed at /? = 5.7(upper) and /3 — 
5.0(lower). These distributions are qualitatively simi- 
lar to the Polyakov- loop distributions evaluated with Dq. 
Such behaviors can be understood intuitively: Di con- 
tains (3fc + 1) wraparound quark loops, and hence corre- 
sponds to the free energy of (S/c-l- 1) light quarks f8|. The 
physical situation is similar to the S/c-quark vacuum (ob- 
tained by evaluating Dq) with one Polyakov loop, where 
(3A:-t- 1) net charges exist. Di^2 can be another candidate 
for an order parameter [8[. 

If we directly generate unquenched gauge configura- 
tions adopting Do instead of D via some adequate algo- 
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FIG. 6: The distributions of Vi computed at /3 = 5.7(upper) 
and (3 = 5.0(lower) are plotted in the complex plain. These 
distributions are qualitatively similar to the Polyakov-loop 
distributions evaluated with Cq. 



in Ref. 



FIG. 4: The distributions of the Polyakov loops evaluated 
with d (upper) and I>2 (lower) ai (3 — 5.0. 



■>-' 



-2 - 




- — ^: D full determinant 



0.1 0.2 0.3 0.4 0.5 



FIG. 5: The phase angles of Do and D obtained as a function 
of chemical potential /i are plotted. They were evaluated with 
one gauge configuration. 



rithms [lil, [la], what is measured is the ratio, 2?i, 2/2^0, 
which will yields 0(1) values. 2?i,2/2'o measured in zero- 
triality sector gives 



(2?„/Po) 



r=o 



/i?[/P„e-^«[^l 



(5) 



which differs in normalization from the original proposal 



/ DUVn 



.-SglU] 



jDUVe-SAU] 



(6) 



The latter can be obtained by measuring 



(2?„/Po 



.T=0 



jDUV^e 



-SAU] 



(2?/2?o) 



T=0 



jDUVe-'^^iu] ' 



(7) 



in zero-triality sector. Pi,2/2?o as well as a Polyakov loop 
evaluated in zero-triality sector also seems usable for the 
clarification of the color confinement. 

If we allow non-zero triality sectors in unquenched 
gauge- updation processes (ordinary updations), the or- 
der parameter should be 2?i^2/2?, which also turns out to 
be the original form proposed in Ref. Q (Eq.([ni)). This 
implementation, however, would not be suitable for such 
purposes, simply because ergodicity is worse. Though 
the configurations that have smaller statistical weights 
De^^o will produce larger I?i.2/2? and Z-^ symmetry in 
(I?i.2/2?) could be finally recovered, such configurations 
would less appear in actual calculations. 

We note here that the inevitable failure of such treat- 
ments in the thermodynamic limit was discussed in 
Ref. 0. Even so, if we directly measure (I?„ )/(!?) in 
quenched backgrounds (like our analyses in this paper 
and the original form in Ref. Q), the ergodicity itself is 
not lost [9[ and this simple prescription could be valid, 
though it is time-consuming. 



C. Spontaneous Zs-symmetry breaking 

The deconfinement phase transition, if it exists, would 
be detected still as the spontaneous Zs-symmetry break- 
ing, as in the quenched case. (Even when the transition 
is crossover, the argument is not invalidated.) We have 
to introduce a small explicit Zs-symmetry breaking ef- 
fect and take the thermodynamic limit before removing 
the explicit breaking, in finite volume systems. Practi- 
cally, it will be workable for the detection of transition 
to evaluate (Rei7'^)3 . Here J7 is a Polyakov loop P or the 
fermionic determinants Pi 2- Especially, the fluctuation 
of 2?i.2 is much suppressed than that of a Polyakov loop 
in zero-triality sector, which can be seen in Fig. [6l[ upper) : 
The values oi'Di^2 are very sharply distributed on the Z3- 
axes, and we can find three "lines" in the complex plain. 
In any case, the use of 2?i,2 seems much more advanta- 
geous than a Polyakov loop, since it gives us more clearer 
signals. 

The reason for the tiny deviations from the Z3 axes 
in Fig. [SJupper) is the discontinuous behavior in D12 
in the deconfined phase. To see this, we compute Pi 
on the gauge configurations where Polyakov-loop's abso- 
lute values are less than 0.4. Such Polyakov-loop's dis- 
tributions are round-shaped both at /3 = 5.7(dcconfined) 
and 5.0(confined). On the other hand, 2?i at /3 = 5.7 
(Fig. El^upper)) still reproduces the three-peak structure, 
which implies 2?i at small \P\ also reflects the phase. 
To have a closer look, we classify configurations into 
three categories, sector 1 (— f < argP < -1), sector 2 
(f < argP < tt), sector 3 (— vr < argP < — f ), in terms 
of the phase angle (argP) of a Polyakov loop. The scat- 
ter plot of Vi in each sector can be found in Fig. [T] 
Though, at /? = 5.7, the distribution of a Polyakov loop 
(|P| < 0.4) is round-shaped, that of Vi is split onto Z3 
axes, which are drawn as three solid lines in the figure. 
This splitting indicates that 2?i_2 "as a function of P" 
quickly changes its value at the boundaries between sec- 
tors in the deconfinement phase whereas they seem to be 
rather smooth functions in the confinement phase. 'D12 
would be discontinuous at the boundaries in the thermo- 
dynamic limit, which is considered as the consequence 
of the spontaneous Za-symmetry breaking. Even if the 
Z3 symmetry in the Zs-symmetrized canonical formula- 
tion is always spontaneously broken in the presence of 
dynamical quarks and the transition is crossover, it will 
be still reflected in 2?i 2 . 



IV. SUMMARY 

We have considered prescriptions that are free from 
the direct charge-screening effects by quark loops and 



enable us to clarify the conflning nature of a vacuum. 
We have tested two candidates for an order parameter, 
a Polyakov loop (P) evaluated in zero-triality sector and 
fermionic determinants (Pi, 2) with non-zero triality. We 
have also individually investigated the distribution of a 

D, evaluated with |P|<0.4 at p=5.7 
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FIG. 7: The distributions of Vi at j3 = 5.7(upper) and 
l3 = 5.0(lower), obtained on the gauge configurations where 
Polyakov-loop's absolute values are less than 0.4, are plotted 
in the complex plain. 



Polyakov loop in each triality sector. Especially, Pi, 2 has 
much smaller fluctuations in comparison with a Polyakov 
loop in zero-triality sector, and seems to well reflect the 
characteristic of a vacuum. Such prescription could single 
out the conflnement nature of a vacuum properly and 
independently of the direct screening effects or the string- 
breaking effects. 
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